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Abstract 

We investigate the magnetic susceptibility of the QCD vacuum with the 1/N C corrections taken 
into account, based on the instanton vacuum. Starting from the instanton liquid model we derive 
the gauged light-quark partition function in the presence of the current quark mass as well as of 
external Abelian vector and tensor fields. We consider the 1/N C meson-loop corrections which are 
shown to contribute to the magnetic susceptibility by around 15 % for the up (and down) quarks. 
We also take into account the tensor terms of the quark-quark interaction from the instanton 
vacuum as well as the finite-width effects, both of which are of order 0(1/N C ). The effects of the 
tensor terms and finite width turn out to be negligibly small. The final results for the up-quarks 
are given as: *0)o — 35 — 40MeV with the quark condensate (iV^V')o- We also discuss the 

pion mass dependence of the magnetic susceptibility in order to give a qualitative guideline for the 
chiral extrapolation of lattice data. 
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I. INTRODUCTION 



The QCD vacuum is known to be one of the most complicated states entangled with 
perturbative and strong nonperturbative fluctuations. In particular, the condensates of the 
quarks and gluons reveal their non-perturbative aspect and character. The quark conden- 
sate is taken as an order parameter associated with spontaneous chiral symmetry breaking 
(S^SB) which is probably one of the most important quantities for low-energy hadronic 
phenomena. Actually, the presence of an external source allows a more profound study of 
the QCD vacuum. For example, the constant electromagnetic field induces another type 
of condensates leading to the nonzero magnetic susceptibility x °f the QCD vacuum defined 
as: 

= e f Xf (i^ f ipf)oF^, (1) 

where ej denotes the quark electric charge with the corresponding flavor /. It is natural to 
have the quark condensate (iipfipf)o in the chiral limit as a normalization factor in the right- 
hand side of Eq.(TjQ), since SxSB is responsible for these quantities, i.e. the quark condensate 
and magnetic susceptibility. The value of Xf was already predicted in the QCD sum rule 
and vector-dominance model 0, 3 B |: Xn(i^Mo = 40 - 70 MeV at the scale of 1 GeV, 
and in the instanton liquid model for the QCD vacuum [5] to ~ 38 MeV and in Ref. [6] 
to 45 ~ 50 MeV. 1 Moreover, Ref. 0] suggested that the magnetic susceptibility x ma y be 
measured in the exclusive photoproduction of hard dijets 7 + N — > (qq) + N. 

In the present work, we want to extend the previous work 0], incorporating the 1/N C cor- 
rections [1, 0] . Since the instanton vacuum explains S%SB naturally via quark zero modes, 
it may provide a good framework to study the x of the light-quark vacuum. Moreover, there 
are only two parameters in this approach, namely the average instanton size p ~ ~ fm and 
average inter-instanton distance R ~ 1 fm. The normalization scale of this approach can 
be defined by the average size of instantons and is approximately equal to p _1 « 0.6 GeV. 
The values of the p and R were estimated many years ago phenomenologically in Ref. ^ 



as well as theoretically in Ref. [111. 1 121. |13| . Furthermore, it was confirmed by various lattice 



simulations of the QCD vacuum 14, 15, 3]. Also lattice calculations of the quark propaga- 
tor 17|, ll8| are in a remarkable agreement with that of Ref. ll|] . A recent lattice simulation 
with the interacting instanton liquid model obtains p ~ 0.32 fm and R ~ 0.76 fm with the fi- 
nite current quark mass m taken into account [19]. Considering the 1/N C corrections, Ref. [9| 



found that the values of the instanton size and average distance should be readjusted as 

p = 0.35fm, i? = 0.86fm (2) 

with F n = 88 MeV and (iip^^o — (255 MeV) 3 required in the chiral limit. The values in 
Eq.Q will be used in the following calculations. 

In order to consider the 1/N C corrections to the magnetic susceptibility of the QCD 
vacuum, we closely follow the formalism developed in Refs. @, 0]. In order to calculate the 
magnetic susceptibility, we first have to find the light-quark partition function Z[V,T,m] in 
the presence of the Abelian external vector and antisymmetric tensor fields. So, we 
first calculate a zero-mode approximated quark propagator S in the instanton ensemble and 



1 In the following we give all numbers for the up-quarks, which are due to assumed isospin invariance 
identical to the numbers for the down quarks. Thus we ignore the index u or /. 
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in the presence of these external fields. Using this propagator, we are able to compute the 
low-frequency part of the quark determinant for the non-zero m with the external Abelian 
vector V and tensor T fields. The relevant techniques have been already developed in 
Refs. @. 

The smallness of the packing parameter 7r(-^) 4 0.1 makes it possible to average the 
determinant over collective coordinates of instantons with fermionic quasiparticles, i.e. con- 
stituent quarks ip introduced. The averaged determinant turns out to be the light-quark 
partition function Z[V, T, m] which is a functional of V and T and can be represented by a 
functional integral over the constituent quark fields with the gauged effective chiral action 
S[^',ip, V,T]. However, it is not trivial to make the action gauge-invariant due to the non- 



locality of the quark-quark interactions generated by instantons. In the previous paper [20 
it was demonstrated how to gauge the nonlocal effective chiral action in the presence of 
the external electromagnetic field and was shown that the low-energy theorem of the axial 



anomaly relevant to the process GG — > 77 is satisfied (see also Refs. [2jJ, |22|). The gauged 
effective chiral action was also successfully applied to various observables of mesons and 
vacuum properties 0, [H, 24. 25l . 



The present work is organized as follows: In Section II, we review the general formalism 
for the magnetic susceptibility of the QCD vacuum. We first show how to derive from the 
instanton vacuum the gauged effective chiral action in the presence of the current quark mass 
as well as of external vector and tensor fields. We then derive the meson propagator for the 
meson-loop corrections to the magnetic susceptibility of the QCD vacuum, considering the 
fluctuation around the saddle-point. In Section III, we discuss the meson-loop corrections 
to the magnetic susceptibility, the contribution of the tensor quark-quark interactions to it, 
and the effects of the finite width of the instanton size, all of which are 1/N C corrections 
to the magnetic susceptibility. In Section IV, we present the final results of the magnetic 
susceptibility of the QCD vacuum and discuss its relevance to the chiral extrapolation in 
lattice QCD, extending the results with larger pion masses employed. In the final Section, 
we summarize the present work and draw conclusions. 

II. GENERAL FORMALISM 

A. Light quark partition function in the presence of external vector and tensor 
fields 

The light quark partition function Z[V, T, m] with the quark mass m and external vector 
Vfj, and tensor T^ v fields is defined as 

Z[V, T,m} = J DA^ e-^ G2 Det(tp + A + im + f + a^T^), (3) 

where is the gluon field and G^ is the gluon field strength tensor. The basic assumption 
of the instanton liquid model is that one can evaluate the integral in a quasi-classical ap- 
proximation, expanding it around the classical vacuum. The first evaluation of the partition 



function Eq. ([3]) was performed in Ref . [28|, |29| in the absence of the external fields and in 



the chiral limit. The main purpose of the present paper is the extension of those results 
to the case of the nonzero quark mass m and external V, T fields, following the method 
given in 0, 0]. For this we split the quark determinant into the low- and high-frequency 
parts according Det = Deti ow Dethi g h and concentrate on the evaluation of Deti ow , which is 
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responsible for the low-energy domain. The high-energy part Dethigh is responsible mainly 
for the perturbative coupling renormalization. 

We start with the zero-mode approximation for the pro pag ator of a quark interacting 
with the z-th instanton. This propagator was considered in 26, 27, 28, 2^]: 

S(= i .LMN (4) 

P + f\i + im p im 
While this zero approximation is good for small values of the current quark mass m, we 



need to extend it beyond the chiral limit as proposed in our previous works |6j, |2Cj, [30|, [31 
as follows: 

Sl = S + So ^^ p'So, (5) 

where 

d = -($oi|^ ^|*oi) = im(<$> m \S p ! \$oi) = im($oi|^ |$oi)- (6) 

The approximation given in Eq.([I)]) allows us to project Si to the correct zero-modes with 
the finite m: 

Si|$oi) = -!-|$oi>, ($0,1^ = ($o,| — . (7) 
im im 

We can write the total quark propagator S in the presence of the whole instanton ensemble 
A and external vector (V) and tensor (T) fields, and the quark propagator with a single 
instanton A, as well as V and T as follows: 

g = \ g = \ fg\ 

p' + A + f + a^T^+im' 1 p 1 + A % +f + cr^T^ + im' KJ 

Here, we assume that the total instanton field A may be approximated as a sum of the single 
instanton fields, A = J2iLi ^ which is justified with the values of p and R in Eq.(j2J). With 
the instanton fields turned off, we write the quark propagator in the presence of the external 
fields and free quark propagator as follows: 

^° = j. i Mr i m — r~ — ' ^° = rr — • (9) 

p + J + a^T^ +im p + im 

We now expand the total quark propagator S with respect to the single instanton field Af. 
S = S + - S ) + J2(Si - S^S^iSj - S ) + ■ ■ ■ . (10) 

i ijtj 

The next step is to expand S with respect to the external fields V and T, and express Si in 
terms of Si. Since we use the zero- mode approximation, the expansion with respect to the 
vector field breaks the gauge invariance. In order to restore it, we introduce the following 
auxiliary field V and gauge connection L,: 

Y^uy+wu-t. (ii) 
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The gauge connection Li is denned as a path-ordered exponent 

Li(x,Zi) = Pexp^iJ^ dy^V^y) 

Li(x,Zi) = j L\(x, Zi)j , (12) 

where z\ denotes an instanton position. The field V/(x, Z{) under flavor rotation if)(x) — > 
U(x)if)(x) is transformed as V-(x,Zi) — > U(zi)V-(x, Zi)U~ 1 (zi). The propagators Si and So 
then have the following form: 

Si = LiSiLi, S> = t + Ji i + y> + a^T» u + im> 

S = LiS> 0i Li, S' 0l = t + K + v^ + im' (13) 
Expanding S' { with respect to Y( and resumming it, we get 

S'i = Si(l + $>V7^D = S' 0i + S' J^y tS'n, (14) 



where 



h = (^(S'^-So^oi), (15) 
(k-h = -($oi|/*S^|$ i>- 



Rearranging Eq. ffTO]) for the total propagator, we obtain 



S = S + SoJ^LT'mo) (Ap + Z^ C Z^ + • • ■) (QoMLfSo 

i,j V / ij 

i,j \ / ij 



where 



Vh = ($od^(^ 1 ^7 1 ¥l$o,)-($o l |^oA rl ^l^), 

T id = (l-<y<$ ^oi^%^o;>, 

T-^ij ^ij^ij = [pi Cij^iji 

Cj = (l-iy)Vy. (17) 
We introduce now the modified zero-mode solution: 

l^o) = -L^o), (18) 
which has the same chiral properties as the zero- mode solution |$o)- Then we get 

s-s = -s J2^}(M (Vpr) l^)(0o,l^ o (19) 
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with 

V + T = fS^. (20) 
The final explicit form for Eq. (|T9|) is written as 

Tr(S -S ) = - X>°.^ ffi)t\M(M |0o,i). (21) 

In order to derive the low- frequency part of the quark determinant, we now introduce a 
matrix operator B(m) defined as follows: 

B(m) tJ = (MySafilfoj) = ($ # (V SoL] 1 ) t\*oj), (22) 
where i,j are indices for the different instantons. Then, we can show that 

k (Det lOT ) = TrlB ( ip + y + T + ~ ) (*» 



= iTr J dm'(S(m') — Sq(i7i')) 
i,j J m 

where Tr denotes the trace over the subspace of the quark zero modes. Thus, we have 

Deti ow [V, T, m] S DetB(m), (24) 

where B is the extension of Lee-Bardeen matrix B [27(| in the presence of the external vector 
and tensor fields V and T. Taking m to be small and switching off the external fields, we 
can show that B turns out to be the same as B to order 0(m). 

Averaging Deti ow over the instanton collective coordinates £, which provides the partition 
function Z^[V, T, m], we can obtain the fermionized representation of the partition function 
in Eq.(|24l: 

Z N [V,T,m] = (Det low [V,T,m\) = (DetB) 

Dil>Dtf exp d^xi)\i> + f + a^ v + im)i/j J] W±ty)\ V>], (25) 



where 



x 



d 4 y (Qlflto&ML-^y^Myj) ■ (26) 

The fermion fields ip ',if> are interpreted as the dynamical quark fields or constituent quark 
fields induced by the zero-modes of the instantons. 

Note that the partition function of Eq. (1251) is invariant under local flavor rotations due 
to the gauge connection L in the interaction term VKt^ While we preserve the gauge 
invariance of the effective chiral action by introducing the gauge connection, we have to 
pay a price: The effective action depends on the path in the gauge connection L. We will 
choose the straight-line path, though there is in general no physical reason why other choices 
should be excluded. However, in Refs. 0, [9| it was shown explicitly that for the magnetic 
susceptibility the path dependence does not come into play. 
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B. Chiral effective action in the presence of external vector and tensor fields 



We are now in a position to derive the relevant partition function for the magnetic 
susceptibility of the QCD vacuum. As exposed in Refs. ^ we first have to average the 
low-frequency part of the quark determinant Deti ow over collective coordinates £±. Here we 
assume a distribution of instanton sizes with vanishing width, i.e. of the form d(p) = S(p—p). 
Having integrated over collective coordinates and having made a exponentiation, we reduce 
the partition function to the following form (for the case Nf = 2): 



r 



Y ± 



P 2 



or 



J, 



a 



fJbV 



J, 



r 



dX + dX^DipDip exp(-r), 



K 

N± In N±+ ^ (ij 

A± 

a 2 det J* + (3 2 det J±, 
f / M 

1 2N 1 



8N r 2N - 1 8N r - 4 



(27) 

(28) 
(29) 

(30) 



'1±75 



if 1 ±75 

T <V — z — W 



If i 

2 ( "/w ± -^ e p,vp\Jp\ 



2 



(31) 

(32) 
(33) 



where the determinant runs over the flavor space tacitly, and K is some inessential constant 
making the argument of the logarithm dimensionless. From Eq. (1301) we immediately see that 
the contribution of the tensor terms is just one of the l/iV" c -corrections. For simplicity, we 
will consider the tensor terms J^ v later. 

Having taken N + = N_ = N/2, having integrated over fermion fields 20, 0, [2^, [3l|, 32 



and having made a bosonization, we end up with the partition function in the presence of 
the external vector and tensor fields: 



Z N [V,T,m} = J d\D<f>exp(-T[V,T,m,\,<l>]) 



(34) 



where 



r[v,T,m,A,$] = Ann — -iv + r^ + r^ 

A 



(35) 



with 



= 2 / d 4 x$ 2 



-Tr In 



P + cr nJT^v + im 



(36) 
(37) 
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Here, $ denotes the meson fields $ = {$ , <&} = {a, 77, cr, 7r}, and $ 2 = $g + $ 2 = d 2 + 
7] 2 + cr 2 + 7r 2 . The P M = p M + eV^, = id^ and Tr stands for the functional trace, i.e. 
jd 4 xtr c tr D tr f . The g is the color factor defined as g 2 = 2N C (N 2 - 1)/(2JV C - 1). The F(p) 
stands for the quark form factor generated by the fermionic zero modes. Though the explicit 
form of F(p) is expressed in terms of the Bessel functions (28j |. we use in the present work 
the dipole type of the form factor for simplicity: 

m - ^r, (») 

where A is the cut-off mass with A ~ 1/p. In fact, the results with Eq. (1381) are very similar 
to those with the original form factor [6j. 

In order to calculate the correlation functions, we follow the general effective action 



approaches [33|, |34|. The effective partition function is now written as 



Z N [V,T,m) =exp(-r eff [V,r,m,A,<$)o]) (39) 

with the effective action T eS [V,T,m, A, ($)o], where vacuum average of the fields ($) must 
be taken solution of the vacuum: 

° r eff [V > r > m,A,($>o]=0 (40) 



and of the coupling A from the saddle-point condition: 

^T eS [V,T,m, A, <$> ] = ^T cS [V,T,m,\, ($) ] = 0. (41) 

C. Vacuum in the presence of external vector and tensor fields 

The instanton vacuum naturally realizes S^SB, which is measured by the quark conden- 
sate (i , 0/'0/}o- The strong quark-quark interaction generated by the instanton vacuum fl26|) 
brings about the non-zero vacuum expectation value of the scalar-isoscalar meson field, i.e. 
a = ($0)0? which can be obtained by minimizing the effective action: 

<9r cff [m, A,a] = dT eS [m,X,a] = 
da ' dX 

Expanding the effective action in one-loop meson approximation, we obtain the following 
expression: 

r cfT [m, A, a] = r [m, A, a] + T™ s {m, A, a], (43) 
where T^ s [m, A, a] represents one-loop meson contribution and is expressed as follows: 



r-KA,a] = iTrln^ 2r ° [m ' A ' $] 



2 V SQiSQj 



(44) 



Here the T is just identical to the effective action in Eq. fl35l) . In Ref. [9j, it is extensively 
studied how to solve Eqs. fj4"2l) and 
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We now find the vacuum average of the meson fields and coupling A in the presence of 
external fields V and T. Though, in fact, their influence appears from the second order, it 
is enough to consider them to the first order in order to calculate the correlation functions. 
Thus, we can safely use a and A derived from Ref . 0] . Then, we can straightforwardly derive 
the meson-loop effective action in the presence of the vector and tensor fields similar to Eq. 



T cS [V, T, m, A, a] = T [V, T, m, A, a] + T™ S [V, T, m, A, a], (45) 

where 



r- ff es [y,T,m,A,a] = ^Spln 



4Sij - -4tt 1 



a 2 P + o^T^ + ifx(P) 
1 



x YiMiP)- — TjM(P) 



(46) 



with n{P) = m + M(P) and M(P) = MF 2 (P). The M denotes the dynamical quark mass 
defined as 

M = ^»^o. (47) 

The Sp stands for the trace over the meson degrees of freedom. 

We have to keep for the self-consistency of the one-meson-loop approximation the meson 
propagators in the leading order (without meson loops) from the equation: 

d I^A = o. (48) 
da 

Then, the meson propagator can be redefined as 

2 V 

n,(,) =* m = ^s(p) + ^(S( P )r,s( P + 9 )r,)- (49) 

III. MESON-LOOP CORRECTIONS TO THE MAGNETIC SUSCEPTIBILITY 
OF THE QCD VACUUM 




FIG. 1: The diagrams which contribute to the magnetic susceptibility in the leading order (first) 
and next-to-leading order (NLO) in 1/N C . The triple line on the left corresponds to the local tensor 
current ipa^ip, the double line on the right represents the vector current ip^^ip. The interaction 
of the vector current with quarks contains both the local and nonlocal terms (see Appendix lAl and 

El). 
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We are now ready for the calculation of the magnetic susceptibility \ from the instanton 
vacuum, taking into account the meson-loop corrections. Diagrammatically, we consider 
three different contributions as shown in Fig. [TJ The first diagram in Fig. [TJ, which is the 
leading-order (LO) contribution in the large N c expansion, has been already calculated in 
Ref. 0: 



x(#V)o = 4iV c 



m 



d A p 

W? W + KpYY (p 2 + m 2 ) 2 



-4N r 



d*p pM{p)F'{p) 
(2tt) 4 (p 2 + /x(p) 2 ) 2 



(50) 



where the quark condensate in the chiral limit plays the role of the normalization. The in- 
frared region of the first integral in Eq. (l50|) brings out the contribution of order 0(m In m) 
which is almost model-independent. In order to calculate the magnetic susceptibility ex- 
plicitly, we first consider two different regions while integrating: < p < 1 GeV and 
1 GeV < p < oo. Then we are able to evaluate the part of the free quark analytically. 
The second integral in Eq. (1501) is related to the nonlocal contribution without which the 
vector current is not conserved, i.e. the Ward-Takahashi identity is broken. It arises from 
the nonlocal quark-quark interaction and is called nonlocal current. 

The meson-loop corrections can be derived by differentiating the mesonic effective action 
given in Eq. fl46p with respect to the external tensor field: 



rrmes 
01 cff 



T=0 



F{P) 



i 



x 



P + i//(P) 



F(P) 



f(p)) r 

/ xy 



a 



[IV 



(51) 



which corresponds to the second and third diagrams in Fig. [TJ Note that here the order 
of the operators should be kept properly, since [P M , P v ] = ieF^ u . Using the Schwinger 
method [H, 36[, we continue to calculate Eq. (l5Tl) . Having carried out a laborious but 
straightforward calculation, we end up with the expression for the meson-loop corrections 
to the magnetic susceptibilities: 



d 4 q 



(52) 



where V* are the vertex functions which are given explicitly in Appendix A and B. The 
Ilj(g) denote the meson propagators given in Eq. p9|) . The final expression for the magnetic 
susceptibility is now written as 



X $ip) = 4iV c 
d 4 p 



d^p n{p)-pMF{p)F'(p) 



AN, 



(2vr) 4 
rn 



(p 2 + /i 2 (p)) 2 



(2tt) 4 (p 2 + m 2 ) 2 (2tt) 4 ' j 
The numerical result for Eq. (l53l) is finally obtained as follows: 



(53) 



Nr. 



0.015 GeV + 5.3 x 10" 
0.415m\ 
GeV J 



m , 
+ ~ — ~ In 



0.007 GeV 



GeV J ' 2tt 2 
0.198m , 

■w 1 1,1 



m y 
GeV/. 



m 
GeV 
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0.038 GeV- (°-^\ - (°-^\ In (- ^ 
+ 0[m\±- 2 ). (54) 



GeV 



IV. CONTRIBUTION OF THE TENSOR TERMS TO THE MAGNETIC SUS- 
CEPTIBILITY OF THE QCD VACUUM 

Since the meson- loop corrections, however, are of order O(l/N ) in the large N c expansion, 
we have to take into account the tensor terms in Eq. (1291) which are also of 1/N C order. It 
is rather tedious to consider the tensor terms [||. In the present Section, we want to show 
briefly how to associate with the tensor terms for the magnetic susceptibility. The details 
can be found in Ref. [9J. 

Introducing the tensor meson fields, the bosonized effective partition function can be 
written as follows: 

Z N [V, T,m] = J dXdXDQDQ^ exp{-T[V, T, m, A, $, VD> ( 55 ) 



where 



T = -N±\n\± + 2[ $2 + -$ 2 



2 



Tr In 



P + T^o^ + im + iV\LF{p) [a^iTi + ^^^a^vA F(p)L~ x 



(56) 



Then, the gap equation and meson fields are changed as follows: 

iVALF(p)(a$ f r< + l^^a^T^F^L- 1 



N 
V~o ' 


2 


*i = 


'-Tr 
4 




1 

-Tr 

4 



P + im + i^XLF{p) (a^Ti + ^^a^Fip)^ 
tV\LF(p)aT t F(p)L- 1 



l\f\LF{p)(3a^T l F{p)L- 1 



4 \p + i m + iV\LF(p)(a^ i T i + f/^^I^F^).^ 1 



(57) 
(58) 
(59) 



In general, it is of great difficulty to evaluate Eqs.(l571)-(|59l. However, if varies very 
slowly, the calculation becomes much simpler. Indeed, for = const, in Eq. f[5"$j) . we get, 
respectively, 



v« c i^ + °(n- (60) 

The constant c\ is proportional to the coupling of the tensor meson to the vector current 
Vfj,, which is just the inverse of the tensor-meson propagator, as will be shown below. Then, 
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we can show that saddle-point values of A and a acquire only (9(V 2 )-order corrections, but 
not 0(V). Substituting these results into the effective action Eq. (l56l) . one can easily get 



dT 



dT, 



fJ,V 



Tr a 



1 



[IV 



T=0 



p + im + iMLF(p)(l + ^^axpT^F^L 



(61) 



with 



a 



N 
ZVr 



+ 0{V 2 ). 



(62) 



It is then straightforward but tedious to calculate the contribution of the tensor terms to 
the magnetic susceptibility: 

mH)o = m f ^ ia^ap + r j^/J^m (63) 



with 



(2tt) 4 (p 2 + fi 2 ( P )y 



2N r 



Cl 



aa J (2tt) 4 (p 2 + // 2 (p)) 5 



rf 4 p MF 2 (p)(^(p)-pMF(p)F'(p)) 
(p 2 + [i 2 (p)) 2 



c \ acrJ J (27r) 4 (p^+M (p)) J 

In Fig. |2l the contribution of the tensor terms is drawn in the dashed box. 



(64) 




FIG. 2: The Feynman diagram for the tensor contribution. 

The numerical results for the contribution of the tensor term are as follows: 

x(^V)o ens ° r = 0.28 x 1CT 3 GeV, ~ .006. (65) 

Thus, the contribution of the tensor terms to the magnetic susceptibility turns out to be 
below 1 %, so that we safely neglect them. 
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V. FINITE- WIDTH EFFECTS ON THE MAGNETIC SUSCEPTIBILITY OF THE 
QCD VACUUM 



We have assumed so far that the width of the instanton-size distribution can be approx- 
imated to zero, i.e.: 

d(p)=5(p-p), (66) 
where all instantons have the same size p. This approximation is justified in the large N c 



limit as shown in Refs. [ll|, [28 

Since we, however, consider the 1/N C corrections to the magnetic susceptibility of the QCD 
vacuum, we also have to take into account the effects of the finite width, which are of order 
O(l/N ) as well. For the numerical evaluation, we take the value for the fluctuation of the 
width as 

V = (p 2 )"(p) 2 -^GeV- 2 , (68) 



which arises from the two- loop size distribution of Ref.[TJJ. In order to take into account the 
effects of the finite width, we first return to Eq. (1261) . in the evaluation of which Eq. (166]) was 
assumed. The integration over p does not change the exponentiation procedure so that we 
still can obtain the standard 2Nf quark-quark interaction in the effective action T. However, 
we should slightly modify the standard procedure to carry out the bosonization. For any 
function J(z, p) we have the identity 

exp J d 4 zdpj 2 (z,p)^j = J D$exp J d 4 zdp<$> 2 + J d 4 zdp$J^J . (69) 

Thus, the effective action is written as 

T = — lnA + 2 / d 4 zdp$ 2 (z,p) (70) 
Vo J 

- Trln (p + a^T^ + im + icVX J d 4 zdpK(z, p)$(z, p)j , (71) 

where c is an inessential constant and K x ^ y (z,p) ~ <j)(x — z,p)(f)(z — y,p). The LO gap 
equations are then expressed as 

N_ = l T[ ( tcV\fd 4 zdp$(p)K(z,p) \ 
V " " 2 1 \^ + im + tc^Jd 4 zdpK(z,p)^(p) J ' 

*o(p) = 7Tr ( ^4 ) , (73 ) 

4 \t + im + icV\J d*zdpK(z,p)%(p) ) 

where $oG°) ^ as ^ e quantum numbers of the a. The unknown variables in Eqs. (I72II73[) are 
the parameter A and the function $(p). In general, these are very complicated equations 
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but can be solved numerically. However, when the width is small, we can expand them 
with respect to 5p 2 . In particular, we expand the functions $(p) and K(p) with respect to 
(p — p), so that we obtain a system of equations for the Taylor coefficients $j and A. Solving 
it, we derive the function $(p) (at least in the vicinity of p). Details of the evaluation can 
be found in Ref. [§]. Finally, we are able to evaluate the corrections of the finite width of 
the instanton size: 

(x(^)o) FW = (-0.00033 + 0.011 (J^L) - 0.096 (^)') • (74) 

The numerical calculation of Eq. (!74|) turns out to be very tiny. The effects of the finite width 
is about 0.6%. Thus, we can safely neglect these effects on the magnetic susceptibility. 



VI. RESULTS AND DISCUSSION 

We have calculated the magnetic susceptibility x(^V)o as a function of the current 
quark mass, with the form factor in Eq. (|38|) employed. For the further discussion it is more 
appropriate to plot the magnetic susceptibility in terms of the square of the pion mass, m 2 . 
Following Ref. j§| both quantities are directly related as 

2 // 1.63\ / K 18.25 13.5577, \ /n/ ' 
m 2 = m ( (3.49 + — - j + m 115.5 + H — — mm) + 0(m 2 ) 

= m (4.04 + 21.587m + 4.52 m In m + 0{m 2 )) [GeV 2 ]. (75) 

Figure [3] shows the magnetic susceptibility V , )o as a function of the square of the pion 
mass, m 2 , with the form factor in Eq. (I38p employed. The long-dashed curve is the result 
of the leading-order contribution in the large N c expansion, which is the same as that in 
Ref. [6J]. The short-dashed curve depicts the meson-loop corrections expressed in Eq. (1521) . 
It is negative for all values of m 2 . Moreover, since the chiral-log term in the leading- 
order contribution is almost cancelled by the meson-loop corrections, the total magnetic 
susceptibility depends on m 2 almost linearly, as shown in Fig. [31 

The meson-loop corrections suppress the magnetic susceptibility by about 15 % near 
the physical value of m 2 = (0.14GeV) 2 = 0.02 GeV 2 . This fact implies that the large N c 
expansion seems to be reliable for them. As a result, the magnetic susceptibility for the up 
and down quarks is: Xu,d = 35 ~ 40 MeV. As for the up and down quarks, the present result 
is comparable to that of Ref. [1, H, HJ . 

Actually, the effect of the meson- loop corrections depends on m 2 . As one can see in 
Fig. [3] at m 2 ~ 0.4 GeV 2 they bring the leading order term down to almost 50%. Thus 
it makes no sense to calculate the susceptibility for larger values of m 2 . Even this limited 
range given in Fig.[3]might be interesting for practitioners of lattice gauge QCD calculations. 
Those calculations suffer generally under the fact that by technical reasons they cannot be 
performed with current quark masses corresponding to the physical pion mass. Usually 
they have to use m 2 > (0.4 ~ 0.6) GeV 2 . In order to obtain a physical observable one has 
to extrapolate from the m 2 used in the lattice calculation to m 2 = (0.14 GeV) 2 . There 
exist different extrapolation techniques as e.g. linear extrapolation, extrapolation by chiral 
perturbation theory finite range extrapolation [42] |43j , extrapolation by 



suitable chiral models etc. Apparently the present calculations suggest in Fig. [3] 

that for the magnetic susceptibility of the vacuum a simple linear extrapolation might be 
sufficient. 
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ml [GeV 2 ] 



FIG. 3: The dependence of the magnetic susceptibility xi^^o on the square of the pion mass, 
m^. The long-dashed curve is a contribution of the leading order result (LO), the short-dashed 
curve is that of the 1/N C -correction (NLO), and the solid curve is the total contribution. Note 
that in the given range the total result depends on almost linearly: The chiral-log terms in the 
leading-order and next-to-leading order results almost cancel each other. 



VII. SUMMARY AND CONCLUSION 

In the present work, we have investigated the magnetic susceptibility of the QCD vacuum, 
based on the low-energy effective QCD partition function from the instanton vacuum. We 
first have constructed the effective partition function in the presence of the external vector 
and tensor fields as well as of the current quark mass. We have made the effective action 
gauge-invariant at the outset. We also have considered the meson-loop corrections which are 
a part of the 1/N C corrections. To be consistent, we have taken into account the effects of the 
finite width of the instanton size distribution as well as the tensor terms of the quark-quark 
interactions, since they are also of order 0(1 /N c ). However, both these effects turn out to 
be negligibly small. 

The meson-loop corrections contribute to the magnetic susceptibility x negatively and 
their chiral-log terms almost cancel those of the leading-order contribution. They contribute 
to the magnetic susceptibility for the up and down quarks by around 15 %. As a result, 
assuming for the up quarks (and identically for the down quarks) m u ^ ~ 5 MeV we obtain 
the magnetic susceptibility x{ii>ty)o = 35 ~ 40 MeV. These results are comparable to the 
estimates of the vector dominance and of the QCD sum rule approaches 0, 0, 0] • It turns 
out that x is almost a linearly decreasing function of the current quark mass m or of m^. 
To know this might be useful for the chiral extrapolation of lattice data for the magnetic 
susceptibility of the QCD vacuum. 
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APPENDIX A: SCALAR VERTICES 



We present the explicit expression for the vertex function for the isoscalar scalar field 
V*(k). The vertex function for the isovector scalar field V*(h) is given as V*(h) = —3V*(k). 
The V*(k) is written as the following expression: 

vm = nJ^- M^ P) F(k +P} [ 12P F {P)F{k+P) ((k +P r + ,(k +P f) 

J (Z7T) 3p(p 2 + M (p) 2 ) 3 ((fc+p) 2 +/i(fc + p) 2 ) 
x (-3p(p)p 2 - p(fc + p)p 2 + 2MF( P ) ( P 2 - fj,(p)fj,(k+p) + k-p) F'(p)p 

+ 2a,(p)V* + p ) - M P )k ■ p) + 2f (p) f ( fc + rt(p 2 + M(p) 2 ) { _ Mk + p)p4 3fc/i(fc + p)p3 

-6/i(fc + p)k ■ PP 2 - 6fcp(fc + p)k -pp + M ((k + p)F{p) (-3p 4 + 2k 2 P 2 - 6p(fc + pfp 2 
-9k ■ PP 2 - 8(k ■ P ) 2 ) F'{p) - P F{k + p) (3p 4 + k 2 P 2 + 6k- PP 2 + 2{k- P ) 2 ) F'{k + p)) 
+/x(p) (-3p 4 - 3fcp 3 - 3fc • PP 2 - 6M(k + p)F{p)ii{k + P )F'( P ) P 2 + 6(fc + p)fj,(k + P f P 
-3kpk ■ p + 4M 2 F{p)F(k + p) (k 2 P 2 - (k ■ P ) 2 ) F'{p)F'{k + p)) 

-3(k + p)p(p) 2 (M (F(p) ( p 2 + k ■ p) F'{p) + Py /{k+p) 2 F(k + p)F'(k + p)) - 2 Pf i(k + p)) ) 
+4 (p 2 + M (p) 2 ) ((fc + p) 2 + ft (k + p) 2 ) (-p (3 M (p) - 2M P F( P )F'( P )) 
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-3 (F(k + p) {pF"(p)(k ■ P ) 2 + (k 2 P 2 + k- PP 2 -(k- p) 2 ) F'{p))) 

- jjr^s (HP) ((k + P )(k 2 + k- P ) 2 F"(k + p) - (fc 4 + 3(fc • p) 2 + (3fc 2 + p 2 ) k ■ P ) F'(k + p)))) 

- (3p/x(p) + 3pp(fc + p) + 2MF(p)k ■ pF'ip)) {F(k+p) {pF'{p) + k ■ P F"( P )) 
- (fc + p)3 ( F W ((k+p)(k 2 + k- p)(p 2 + k ■ p)F"(k+p) 

-(p 4 + 2k 2 p 2 + (k-pf+(k 2 + 3p 2 )k-p)F'{k+p)))))] . (Al) 



APPENDIX B: PSEUDOSCALAR VERTICES 



We present the explicit expression for the vertex function of the pseudoscalar isoscalar 
field V*(k). The vertex function of the pseudoscalar isovector field V*(k) is given as V*(k) = 
—3V*(k). The V*{k) is written as the following expression: 



V^k) = N c jp- MF( P )F(k +P ) U 2p F {p )F(k + P )((k +P ) 2 + ,(k +P ) 2 ) 

J (Z7T) 3p(p 2 + M (p) 2 ) 3 ((fc+p) 2 +p(fc + p) 2 ) 
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(-3[i(p)p 2 + n(k + p)p 2 + 2MF(p) (p 2 + fj,{p)fj,(k + p) + k-p) F'(p)p - 2n(p) 2 fi(k + p) - 3fi{p)k ■ p) 

+ 2F(p)F(k+p) (p +/x(p) ) + + ^ + + 6fi(k+p)k ■ pp 2 + 6kn(k + p)k ■ pp 

(k+p) 

+M ((k + p)F{p) (-3p 4 + 2k V - 6^(k + pfp 2 ~ 9k ■ pp 2 - 8(k ■ p) 2 ) F'(p) 

+ pF(k+p) (3p 4 + k 2 p 2 + 6k ■ pp 2 + 2(k ■ p) 2 ) F'(k+p)) - fi(p) (3/ + 3kp 3 + 3k ■ pp 2 

-6M (k + p)F(p)fj,(k + p)F'(p)p 2 - 6(k + p)^{k + pfp + 3kk ■ pp 

+AM 2 F(p)F(k + p) (k 2 p 2 - (k ■ p) 2 ) F'{p)F'{k + p)) 

+3{k + p)v(p) 2 (M (p(k + p)F(k + p)F'{k + p) - F{p) (p 2 + k- p) F'(p)) - 2 Pf i(k + p))) 
+4(p 2 + M (p) 2 ) ((fc+p) 2 +/i(fc+p) 2 ) 

p(3/i(p) - 2MpF(p)F'(p)) (F(k+p) ( P F"(p)(k-p) 2 + (k 2 p 2 + k ■ pp 2 - (k ■ p) 2 ) F'(p))) 

Tj^i (P(P) ((* + p) (k 2 + k ■ p) 2 F"{k +p)- (k 4 + 3(k ■ p) 2 + (3k 2 +p 2 )k- p) F'{k + p)))) 
(3pfi(p) - 3pn(k +p) + 2MF(p)k ■ P F'(p)) (F(k + p) ( P F'(p) + k ■ P F"(p)) 
(fc * 3 (F(p) (V(k+P) 2 (k 2 + k-p)(p 2 + k- p)F"{k + p) 

~(p 4 + 2k 2 p 2 + (k-p) 2 + (k 2 +3p 2 )k-p)F'(k + p)))))] . (Bl) 
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